INSTANTANEOUS MEASUREMENT OF SIGNAL POLARIZATION 
BACKGROUND OF THE INVENTION 

[0001] The present invention relates to apparatus and methods 
for determining the polarization of electromagnetic signals. 
[0002] Electromagnetic signals such as radio waves and light 
have a property referred to as polarization. Radar operates by- 
transmitting an electromagnetic signal to a target and comparing 
the signal reflected from the target with the transmitted signal. 
In modern electronic warfare, targets avoid detection from enemy 
radar by using various countermeasures such as, jamming an enemy 
radar signal impinging on the target with a signal denying range 
information to the enemy and creating false reflected signals to 
deceive the enemy radar system. To be effective, the signals 
created by the countermeasure system should have characteristics 
such as polarization corresponding to the signal characteristics 
expected by the enemy system as, for example, characteristics of 
the return signals expected by an enemy radar system. In some 
cases, the enemy radar may change its signal polarization 
rapidly. Such a radar system is referred to as a "polarization 
agile." If the enemy radar is polarization agile, the 
countermeasure system must be capable of determining the 
polarization of the transmitted signal rapidly, so that the 
countermeasure system can change the signals which it emits. For 
example, a jamming system carried on an aircraft and intended to 
defeat a polarization agile enemy radar system should determine 
the polarization of the incoming radar signal and alter the 
polarization of the jamming signal accordingly. If the jamming 
system does not do this, the jamming signal will not match the 
polarization of the return signals from the aircraft. The enemy 
radar receiver can reject the jamming signals and acquire 
meaningful return signals . Delay in measuring the incoming 
signal polarization can allow the enemy system to acquire 
meaningful return signals for a sufficient time to find the 
position of the aircraft. Conversely, where a radar or 
communications system must overcome enemy jamming, it is 
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desirable to measure the polarization of the jamming signal and 
transmit the radar or communications signal with a different 
polarization . 

[0003] However, traditional polarization measuring techniques 
do not provide polarization measurements rapidly enough to 
counteract a polarization agile enemy system. Just as the 
receiving system becomes accustomed to one polarization, the 
enemy system changes polarization. 

[0004] At a given point in space along the path of an 
electromagnetic wave and at a given instant in time, an electric 
field points in a particular direction, denoted by a vector, E . 
This vector is perpendicular to the direction of travel of the 
signal or "propagation vector." The polarization of an 

electromagnetic wave is described by the orientation of the 
electric field vector and the manner in which this vector varies 
with time. 

[0005] The polarization vector can be split into components E x 
and E y along orthogonal x and y axes perpendicular to the 
direction of travel of the electromagnetic wave. The component 
along the x axis commonly is referred to as the "horizontal" 
component, whereas the component along the y axis is referred to 
as the "vertical" component. Although these terms are used 
herein, it should be appreciated that these directions may be 
arbitrary directions unrelated to the normal gravitational frame 
of reference. At any given point in space, E x and E y vary with 
time. For example, for a sinusoidal wave having frequency ©, 
E x =Asin(cot) and E y =Bsin ( (cot ) +a) , where t is time, a is a phase 
difference and A and B are the magnitudes of the E x and E y 
components. When the E x and E y components are in phase (a =0), 
the electric field is linearly polarized. In this condition, the 
electric field vector at a given point always lies on the same 
plane. When the E x and E y components are out of phase (a ^0) , 
elliptical polarization results. When the E x and E y components 
of an elliptically polarized electromagnetic signal are of equal 
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magnitude (A=B) and are 90° or 27 0° out of phase, the signal is 
said to be circularly polarized. 

[0006] To measure signal polarization, a dual-aperture 
(polarized) antenna and a device known as a polarimeter are 
required. The dual-aperture antenna provides one electrical 
signal V h representing the E x or horizontal component of the 
electric field of a signal impinging on the antenna, and another 
electrical signal V v representing the E y or vertical component of 
the electric field of the same signal. These signals typically 
are amplified and filtered separately in a dual-channel receiver 
before passing to the polarimeter. The polarimeter compares 
these signals to determine their relative magnitudes and the 
phase difference between them. 

[0007] A prior art analog polarimeter is shown in FIG. 1. The 
horizontal signal V h is supplied to one input of a four port 
directional coupler 2 00 of a type referred to as a "hybrid." The 
vertical signal V v is supplied to the input of a phase shifter 
202 which applies a known phase shift $ to that signal. The 
phase-shifted signal is supplied to another input of the 
directional coupler 200. The coupler 200 provides a signal at a 
first output 204 representing the coupled power output or sum of 
the input signals supplied to the circuit, and also provides a 
signal representing a specific phase shift between the input 
signals at a second output 206. In this prior art example, 180°. 
The first or sum output 204 of circuit 200 is supplied to the 
input of a further phase shifter 208 which applies a known phase 
shift y- The output of this phase shifter is connected to one 
input of another directional coupler 210, which is similar to the 
first 200. The second or difference output 206 of combining 
circuit 2 00 is connected directly to the other input of combining 
circuit 210. Thus, when time-varying V v and V h signals are 
applied to the polarimeter, one time-varying output signal, 
referred to as the A signal appears at the difference output 212 
of coupler 210. Another time-varying output signal referred to as 
the £ signal, appears at the sum output 214 of coupler 210. The 
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output signals are supplied to a dual-channel receiver and 
logarithmic amplifier 216 which monitors the amplitudes of these 
signals and provides a signal representing a ratio between their 
amplitudes. This ratio signal is supplied to a null adaptive 
tracker 218, which adjusts the phase differences cj> and y applied 
by the phase shifters to achieve a null condition as discussed 
below . 

[0008] The relationships between the A and S output signals 
appearing at outputs 212 and 214 and the input signals V v and 
are referred to as the "transfer functions" of the polarimeter. 
These transfer functions depend on the phase shifts <)> and y 
applied by phase shifters 202 and 208. Conversely, there is a 
relationship between the transfer functions which yield output 
signals with particular characteristics and the phases and 
amplitudes of V v and V h . Stated another way, there is a 
relationship between the phase shifts and y which yield 

particular output signal characteristics and the phases and 
amplitudes of the input signals V v and V h , 

[0009] In particular, for the components illustrated in 
|A| 

Fig. 1, the ratio t-j between the amplitude |Ajof the A 
H 

output signal and the amplitude |S|of the S output signal will be 
at a minimum or null condition when: 

y = 2taiT'p-j , and (1) 
Where : 

a is the amplitude of the horizontal component V h ; 
b is the amplitude of the vertical component V v ; and 
a is the phase difference between these components. 

[0010] Solving for the amplitude ratio — and phase difference 
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tanjj-j, and (3) 



(4) 



A, 371 
a = 0 — — . 

[0011] Thus, the parameters which characterize the 

polarization of the signal, such as the amplitude ratio ^ and 

phase difference a between the components of the signal can be 
found from the phase shifts <j> and y which yield the null 

condition or minimum ratio j^l . Tilt angle, t, of an 

n 

elliptically polarized signal is also derivable from the 
polarimeter phase shift angles y and (j) at the null condition as 

t = ^-tan" 1 [tan (2y) cos ((() -^-)]. ( 5 ) 

[0012] In operation, tracker 218 sets phase shifter 208 to 
hold y constant at an arbitrary value and adjusts phase shifter 
202 to vary <|> in an iterative or trial-and-error process 

until the output ratio j^j is at a minimum for the arbitrary 

value of y- Tne tracker 218 then holds (|> constant and 
adjusts phase shifter 208 to vary y in a further iterative 
process until the true minimum or null condition is found. 
[0013] Other known analog polarimeters use different networks, 
typically including phase shifters and mixers. However, the 
overall principle of operation is the same. The transfer 
function or functions of the polarimeter is adjusted iteratively 
to yield output signals having predetermined characteristics, and 
the polarization of the signal is found from the transfer 
function or functions which yield those characteristics. 
Polarimeters of this type can provide accurate measurements of 
signal polarization. However, they require considerable time to 
perform the required iterations. 
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SUMMARY OF THE INVENTION 

[0014] One aspect of the present invention provides apparatus 
for determining the polarization of a signal having vertical and 
horizontal components. Apparatus according to this aspect of the 
invention includes a first register for storing a first series of 
discrete values representing a series of samples of the 
horizontal component and a second register for storing a second 
series of discrete values representing a series of samples of the 
vertical component. The different values in each series can be 
denoted by an integer index. For example, a series of N 
horizontal sample values can be represented as A ( 1) , A (2 ) ...A (N) or, 
generically, A() where the parenthetical expression represents 
the index, and the vertical samples values can be represented by 
B() . 

[0015] The apparatus also includes a plurality of sets of 
calculation elements which are arranged to combine values in the 
first and second series with one another so as to produce one or 
more output values. Typically, the calculation elements are 
arranged to operate cyclically, so that a reference value in one 
series of sample values is combined with other values in the 
other series, in the same series, or both during each cycle of 
such set, to yield one or more output values for each cycle. 
Operation of each set of calculation elements through numerous 
cycles, using different values in one series of input values as 
the reference value, yields one or more series of output values. 
Each set of calculation elements has one or more transfer 
functions specifying the manner in which different samples in the 
series are combined with one another. These transfer functions 
typically include one or more integer offsets specifying the 
differences between the index of the reference value used on a 
particular cycle and the index of each other value to be combined 
with the reference value on that cycle. Preferably, the transfer 
functions used by different sets of calculation elements are 
different from one another. At least some of these different 



sets of calculation elements are arranged to operate in parallel 
with one another. 

[0016] The apparatus also includes an evaluation circuit 
connected to the various sets of computation elements. The 
evaluation circuit is arranged to compare one or more 
characteristics of the series of output values generated by the 
various sets of calculation elements with one or more preselected 
characteristics, and to select the series having characteristics 
corresponding to preselected characteristics. This selection 
inherently identifies the set of computation elements which 
provided such series, and thus identifies the offsets used in the 
transfer function of that set. The identified offsets provide 
information about the polarization of the signal. 

[0017] Operation of an individual set of calculation elements, 
with particular offsets, is analogous to operation of an analog 
polarimeter with particular phase delays. However, because 
numerous sets of calculation elements operate in parallel with 
one another, the need for iteration is reduced or eliminated. 
Only a small number of cycles are required to determine signal 
polarization. Polarization measurement can be accomplished 
rapidly, even where many different offsets are used to provide a 
fine phase resolution as required for a high-accuracy 
polarization measurement, and even where each series of sample 
values includes thousands of sample values. 

[0018] For example, a polarimeter according to one embodiment 
of the present invention includes numerous sets of calculation 
elements in the form of adders connected to the registers to 
provide transfer functions of the form: 

A (k, i, j) = A (k) - B(k + i) - [A(k + j) + B(k + i + j ) ] , (6) 
s (k, i, j) = A (k) - B(k + i) + [A(k + j ) + B(k + i + j)], (7) 
In these functions, k is the index of • the reference sample A(k) 
whereas i and j are the offsets and A() and S() are two series of 
output values produced by each set of calculation elements. 
These particular transfer functions are analogous to the transfer 
functions of the analog polarimeter discussed above. For series 
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of input sample values A() and B ( ) corresponding to particular 
input signals V v and V h , a set of these transfer functions 
with particular offsets i and j will provide series of 
output values A() and S() corresponding to the A and S output 
signals of the analog polarimeter of Fig. 1 using particular 
phase shifts (j) and y. Differences in j and i have effects 
analogous to differences in $ and y, respectively. However, 
because the various sets of calculation elements operate in 
parallel with one another, it is unnecessary to vary each phase 
shift iteratively. 

[0019] For example, in an ideal polarimeter according to this 
embodiment of the invention, a first group of N sets of 
calculation elements having transfer functions with different 
values of j but with the same value of i are operated in 
parallel, and the series of output values which yields the lowest 

amplitude ratio is selected to thereby select a value of j. 

That selected value of j is used in operation of a second group 
of M sets of calculation elements each having transfer functions 
with the selected value of j but with different values of i. 
Here again, the series of output values which yields the lowest 

W i s selected, which in this case results in 

selection of a value i. The selected values of i and j 
can be converted into the polarization parameters of the 
signal, such as the as the amplitude ratio ^ between the 

vertical and horizontal components , the phase difference a and 
the tilt angle i. Where each series of input sample values A() 
and B() includes K samples, the polarization of the signal can be 
determined in slightly more than 2K clock cycles. 

[0020] Most preferably, each set of calculation elements is 
associated with one or more characteristic-calculation circuits, 
and the characteristic-calculation circuit associated with each 
set of calculation elements operates in parallel with the 



calculation elements of that set, and in parallel with the 
characteristic-calculation circuits associated with other sets. 
For example, where the characteristics of the output series which 
are examined include amplitude, the characteristic-calculation 
circuits associated with each set of calculation elements may- 
include one or more accumulators each arranged to add an output 
value calculated on each cycle of the calculation elements to a 
total. For example, where each set of calculation elements 
provides a A() value and a S() output value on each cycle, the 
characteristic-calculation circuits associated with each set may 
include one accumulator for adding the A ( ) value produced on each 
cycle to a total £A and another accumulator for adding the S() 
value produced on each cycle to a total SS. Essentially, no 
additional time is required to calculate the amplitudes of the 
output signals from the various sets of calculation elements. 
[0021] In a particularly preferred arrangement, the registers 
used to hold the series of input sample values include shift 
registers. The different offsets used in the transfer functions 
of the various sets of calculation elements are established by 
connecting the different calculation elements to different taps 
of the shift registers. The reference value used in each cycle 
of the calculation elements is changed by clocking the data 
through the shift registers. 

[0022] A further aspect of the present invention provides 
methods of determining the polarization of a signal from a series 
of horizontal input sample values and a series of vertical input 
sample values. Methods according to this aspect of the invention 
desirably start with a calculating a plurality of series of 
output values, using a plurality of sets of transfer functions. 
Most preferably, calculations using at least some of the sets of 
transfer functions are conducted in parallel with other 
calculations using other transfer functions. As discussed above 
in connection with the apparatus, the transfer functions 
represent combination of samples in the two series with a 
reference sample value. Each transfer function includes one or 

9 



more integer offsets specifying the differences between the index 
of the reference value used on a particular cycle and the index 
of each other value to be combined with the reference value on 
that cycle. The transfer functions in different sets desirably 
include different offsets. The method desirably further includes 
the step of evaluating one or more characteristics of the series 
of output values computed using the transfer functions of the 
different sets according to a predetermined criterion, and 
selecting the series produced by one sets of calculations based 
on such evaluation. As discussed above in connection with the 
apparatus, this selection implicitly selects one set of 
calculations and thus selects one set of offsets, from which the 
polarization characteristics of the signal can be determined, 
calculation of numerous series of output values in parallel with 
one another minimizes the need for iteration. 

[0023] Particularly preferred methods according to this aspect 
of the invention, one or more characteristics of each series of 
output values are calculated in parallel with calculation of the 
output values themselves. Methods according to this aspect of 
the invention can provide rapid polarization measurements. 
[0024] Other objects and advantages of the system and method 

will become apparent to those skilled in the art after reading 
the detailed description of the preferred embodiment. 
BRIEF DESCRIPTION OF THE DRAWINGS 

[0025] FIG. 1 is a system block diagram of a prior art 
polarimeter . 

[0026] FIG. 2 is a functional block diagram of apparatus in 
accordance with one embodiment of the present invention. 
[0027] FIG. 3 is a more detailed block diagram of a portion of 
the apparatus shown in FIG. 2. FIG. 3 is presented as two parts, 
FIGs. 3a and 3b. 

DETAILED DESCRIPTION OF THE INVENTION 

[0028] The embodiments will be described with reference to the 
drawing figures where like numbers represent like elements 
throughout . 
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[0029] Shown in FIG. 2 is an apparatus in accordance with one 
embodiment of the present invention is arranged to determine the 
polarization of a received radar pulse. The apparatus includes a 
sample preparation circuit 14 which recovers series of discrete 
digital sample values representing the vertical and horizontal 
components of the signal, and the polarimeter 17 which operates 
on these series to determine the polarization of the signal. 
[0030] The sample preparation circuit includes a dual-aperture 
antenna 16 having two outputs 19h and 19v for the horizontal and 
vertical signal components. A dual channel receiver 21 

separately downconverts the horizontal 1% and vertical signal 
component 19 v to a predetermined intermediate frequency. Because 
the enemy radar signal may change frequency rapidly, receiver 21 
desirably is a frequency-agile receiver controlled by 
conventional frequency-detection circuitry (not shown) . The 
intermediate-frequency horizontal and vertical component signals 
are supplied by receiver 21 to low pass filters 23 h and 23 v , 
respectively. The filtered components are separately digitized 
by analog-to-digital (ADC) converters 25 h and 25 v . 

[0031] The sampling rate of the ADC's 25 h , 25 v , the 
intermediate frequency used by receiver 21 and the pass-band of 
the low pass filters 23 h , 23 v are selected so that the ADCs 25 h , 
25 v will be able to digitize the highest expected frequency 
output from the receiver 21 and filters 23 at a sampling rate 
above the Nyquist rate for that frequency, and desirably at a 
rate of 20 samples/cycle. The ADCs 25 h , 25 v desirably quantize 
the filtered IF signal components with a minimum resolution 
of 8 bits yielding 256 quantization levels. For operation with 
incoming radar signals at about 6 to about 18 GHz, the 
intermediate frequency may be on the order of 10 MHz, and hence 
each ADC 25h,25v operates at a sampling rate on the order of 200 
MHz, i.e., 200 million 8-bit digital sample values per second. 
The ADC's may be converters of the type disclosed in commonly- 
assigned United States Provisional Patent Application Serial 
Number 60/164,947, filed on November 12, 1999, the disclosure of 
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which is incorporated by reference herein. A discussion of the 
conversion system and method is beyond the scope of this 
disclosure . 

[0032] The sequences of vertical and horizontal component 

sample values produced by ADCs 2 5h and 2 5v are stored in received 
signal component memories 27h and 27v. These sequences are 
expanded or "stretched" by oversamplers 2 9h and 29v, which may 
use a conventional interpolation process to insert numerous 
intermediate sample values between each pair of actual sample 
values in each such sequence. This yields a series of horizontal 
sample values A() and a series of vertical sample values B(), 
which are stored in oversampler memories 31h and 31v, 
respectively. In this notation, the parenthetical expression 
denotes an integer index, i.e., the first value in the horizontal 
series is denoted as A(l), the second sample is denoted as A(2), 
and so on. 

[0033] The quality of the polarization measurement produced by 

polarimeter 17 depends upon the phase resolution of the system, 
which in turn depends upon the number of samples per cycle of the 
original signal present in each series. The quality of the 
polarization measurement is referred to as the null 
depth, expressed in decibels (dB) . The relationship 

between phase resolution and null depth is 



null (in dB) 



where 8\\r is phase resolution in degrees. Typically, a 
null depth on the order of 40 dB is desired, which 
implies a phase resolution of 0.5 degrees or better. To 
achieve this phase resolution, each series of sample 
values A() and B() should include 720 samples per cycle. 
Stated another way, within each series A() and B(), each 
increment in the index corresponds to a given phase 
delay, d<j>, which is 0.5 degrees in the case where each 
series includes 720 samples per cycle. Sample values 
with the same index in the two series A() and B() 
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desirably represent portions of the original horizontal 
and vertical component signals occurring at the same 
time, i.e., A(l) and B(l) represent the horizontal and 
vertical components of the signal at the same moment, 
A(2) and B(2) represent the horizontal and vertical 
components at the next moment, and so on. 

[0034] Each series of sample values should include data 
representing more than one cycle of the original signal. In a 
preferred embodiment, each series of sample values represents 
eight full cycles of the original signal, and hence includes 
5,760 (8 times 720) individual sample values. The polarimeter 17 
operates on these series of sample values. 

[0035] The polarimeter includes a first group 37 of N x 
"channels" or sets of computation elements. As further discussed 
below, each set operates with samples having different offsets or 
differences in index. Certain offsets used in the various sets 
differ from one another. For maximum phase resolution, the 
difference in offsets between two sets should be one index value. 
The differences in offsets among all of the sets used in this 
first group should total at least one full cycle or 3 60°. This 
arrangement is used in the embodiment illustrated. A difference 
in offset of one index value corresponds to 0.5°, and hence 720 
sets are required to span the full 360° range. Thus, in this 
embodiment Ni is 720. The functional interrelationships of the 
components constituting this group " of channels and associated 
elements are shown in Fig. 3a. A first-stage horizontal sample 
shift register 102 and two first-stage vertical sample shift 
registers 104 and 106. Each shift register is a conventional 
device defining Ni memory locations arranged in sequence from an 
input or upstream end to an output or downstream end. Each shift 
register has an output or tap associated with each memory 
location. In the conventional manner, each shift register is 
arranged to operate cyclically. In each cycle of operation, the 
value at each memory location is supplied through the output 
associated with that memory location and shifted to the next 
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downstream location. Also, on each cycle of operation, a new 
value is accepted or "clocked into" the first memory location, 
whereas the value at the downstream-end memory location is 
shifted out of the register. Thus, each register can accept a 
sequence of sample values and deliver different sample values in 
the series from the various outputs or taps. For example, 
register 102 will receive the A() series of horizontal sample 
values from memory 31h. On any given cycle of operation, the 
index of the value A() delivered from the most upstream output 
will be lower than the index of the value A() delivered from an 
output at position downstream along the register. The difference 
in index corresponds to difference in position along the 
register. Difference in index also corresponds to a difference 
or delay in phase. 

[0036] Register 104 will receive the series of vertical sample 
values B() from memory 31v. The input of register 106 is 
connected to an output of register 104 i positions downstream 
from the input of register 104. Thus, on any given cycle the 
index of the value B() in each position of register 106 will be 
offset by i from the value in the corresponding position of 
register 104. The value of i used is entirely arbitrary; it may 
be any integer between 1 and Ni. All of the registers operate in 
synchronism, so that when a particular value A(k) is in the 
reference position at the upstream end of register 102 when the 
corresponding or same-index value B(k) is at the input or 
upstream end of register 104. 

[0037] As mentioned above, the first group of computation 

elements includes Nl sets of computation elements. A typical set 
108j of computation elements includes a first adder 110 having a 
positive input connected to the reference output at the input or 
upstream end of horizontal sample register 102 and also having a 
negative input connected to the output of the vertical sample 
register 104 i positions downstream from the input end of that 
register. Thus, in a cycle when sample A(k) is at the reference 
value position, at the upstream end of register 102, adder 110 



will receive sample A(k) through the positive input and sample 
B(k+i) through the negative input, and will provide an output 
equal to A(k)+B(k+i) . 

[0038] The same set 108j also includes a second adder 112 
having a positive input connected to an output of horizontal 
sample register j positions downstream from the reference 
position or input end of that register and having another 
positive input connected to an output of register 106 j positions 
downstream from the input or upstream end of that register. 
Thus, on a cycle when sample A(k) is at the reference value 
position at the upstream end of register 102, the second adder 
114 112 will receive samples A(k+j) and B( k+i+j) will produce an 
output equal to [A (k+j ) +B (k+i+j ) ] . 

[0039] Set 108j also includes a third adder 114 having a 
positive input connected to the output of the first adder 110 and 
a negative input connected to the output of second adder 112. 
Thus, on a cycle when sample A(k) is at the reference value 
position at the upstream end of register 102, the third adder 114 
will yield an output sample value A(k,i,j), where A(k,i,j)= A(k)- 
B(k+i)- [A(k+j)+B(k+i+j) ] 

[0040] Set 108j also includes a fourth adder 116 having 
positive inputs connected to the outputs of the first adder 110 
and second adder 112. Again where value A(k) is at the reference 
value position, the fourth adder 116 will yield an output sample 
value E(k,i,j) where Z(k,i,j)= A(k)-B(k+i)+ [A (k+j ) +B (k+i+j ) ] . 
[0041] Thus, as successive horizontal sample values A(k) and 
vertical sample values B(k) are clocked through the registers, 
the first set of adders will produce two series of output sample 
values S(k,i,j) and A(k,i,j). 

[0042] Every other set of computation elements 108 is 
identical to set 108j discussed above, except that the value of 
the offset j used for each set is different. For example, for 
set 108 N1 ,the value of j is equal to N 2 and hence the 
inputs of second adder 112 are connected to outputs of 
registers 102 and 106 Ni positions from the upstream ends 
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of these registers. Accordingly, as successive values 
A(k) and B(k) are clocked through the registers, each set 
108i through 108 N i will produce a two series of output 
values E(k,i,j) and A(k,i,j) as discussed above, using the 
same value of offset i but different values of offset j . 
[0043] A pair of accumulators 120 and 122 is associated with 
each set of adders 108. For example, the first accumulator 120j 
associated with set 108j receives the series of output sample 
values A(k,i,j) from the third adder 114 of that set, and adds 
the absolute value |A(k,i, j) | of the sample value for each cycle 
to a total SAj . Thus, the total accumulates over successive 
cycles, i.e., over the various values of k. Likewise, the second 
accumulator 122j associated with set 108j receives the series of 
output sample values E(k,i,j) from the fourth adder 116 of set 
108 j, and accumulates the total ESj equal to the sum of the 
absolute values |S(k,i,j)|. The accumulators associated with the 
various sets are identical, but accumulate totals for different 
series of output sample values resulting from different values 
of j. 

[0044] The accumulators 120 j, 122 j associated with set 108 j 
are connected to a pair of logarithm-calculating circuits 124 j, 
12 6 j respectively. When actuated by a control signal, the 
logarithm-calculating circuits associated with will calculate the 
logarithms log(SAj) and log (ESj) . Similar logarithm- 

calculating circuits 124 and 126 are connected to the 
accumulators associated with the other sets 108. The 
logarithm-calculating circuits associated with set 1 0 8 j 
are connected to a difference and multiplication circuit 
128j arranged to calculate 2 0 [ log (SAj ) -log (SSj ) ] , or 
2 Olog ( SAj /SS j ) . A similar difference and multiplication 
circuit 128 is associated with each set 108. 

[0045] In operation, after a series of input samples 
A() and B() have been clocked through the registers, the 
logarithm circuits 124, 126 and the difference and 
multiplication circuit 128 associated with the various 
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sets are actuated. The difference and multiplication 
circuits yield Ni values forming an Ni component vector 
130 of the form 2 Olog ( 2 A1/SE1 ) ... 201og (SAj/SSj ) ... 
201og (SA n i/SSni) . 

[0046] A comparator circuit 132 is connected to the 
outputs of the difference and multiplication circuits 
128, so that it receives vector 130. Comparator circuit 
132 compares the components of the vector with one 
another to find the smallest component. That component 
represents the smallest ratio (EAj/ESj). The comparator 
thus selects the value of j associated with that 
component and outputs that value as the selected value J. 
[0047] The horizontal and vertical sample values A() and B() 
clocked out of shift registers 102 and 104 are routed to buffers 
300 and 301 (Fig. 3b), respectively and stored for use in a 
second group 41 of parallel processing channels or sets. 
Horizontal sample buffer 300 is connected to the input of 
a second-stage horizontal shift register 302, and the 
vertical sample buffer 301 is connected to the input of a 
first second-stage vertical shift register 304. These 
registers are identical to the corresponding registers 
102 and 104 used in the first stage. A switching network 
305 is arranged to make an input connection to any one of 
the outputs of vertical shift register 304. The 
switching network has an output connected to the input of 
a second vertical shift register 306. Switching network 
305 is responsive to the J output of comparator 132, and 
receives the value of J selected by the comparator. The 
switching network makes its input connection to the 
output of register 304 J positions downstream from the 
input or reference end of the register. Thus, as 

successive vertical sample values are clocked through the 
register 304, a series of vertical sample values offset 
by J values will be clocked through register 306. Stated 
another way, when value B(k) appears at the reference 
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position at the upstream end of register 304, value 
B(k + J) will appear at the upstream end of register 306. 
A similar switching network 307 is connected to register 
302. This switching network is also responsive to the J 
output of comparator 132 in the first stage, and makes an 
input connection J positions downstream from the upstream 
or reference position of register 302. Thus, when 

horizontal sample value A(k) appears at the reference 
position of register 302, sample value A(k+J) appears at 
the output of switching circuit 307 . 

[0048] The second group of calculation elements 41 
includes N 2 sets 308 of calculation elements operate. 
Here again, different sets use different offsets. 
However, in this second stage, the differences in offsets 
need only span a range of 180°. Preferably, the same 
one-index difference in offsets between sets is employed. 
Here again, each one-index difference in offsets 
corresponds to 0.5°. Thus, 360 sets 308 are used in this 
stage, i.e., N 2 =360. Except as discussed below, each set 
of calculation elements 308 is identical to a set of 
calculation elements 108 discussed above with reference 
to Fig 3a. However, the first adder 310 of each such set 
has its negative input connected to a different position 
along vertical sample register 304. For example, the 
negative input of first adder 310i in set 308i has its 
negative input connected to an output of register 304 i 
positions downstream from the input end of the register. 
Here again, the first adder 310 of each set has its 
positive input connected to the reference position or 
upstream end of horizontal sample register 302. Thus, 
the first adder 310 of each set will yield A(k)-B(k + i), 
and but using different values of i in each set. The 
second adder 312 of each set 308 has one positive input 
connected to the output of switching network 307 and 
another positive input connected to an output of register 

18 



I 



306. Different outputs of register 306 are used for the 
second adders of different sets. An output i positions 
downstream along register 306 is connected to the second 
adder of the ith set 308i. Thus, the second adder 312 of 
the ith set will yield A { k+J) +B ( k+i+J) . The third adder 
314 and fourth adder 316 of each set are identical to the 
third and fourth adders of sets 108. Thus, on successive 
cycles, the third adder 314 of each set will yield a 
series of values A {k, i, j ) computed as discussed above 
except that for all sets 308 j is the same and is equal 
to J and i is different for the different sets. 
Likewise, the fourth adder 316 of each set 308 will yield 
a series of values S(k,i,j) using different values of i 
for each set but using the same value of j (j=J) for all 
sets . 

[0049] A pair of accumulators 320 and 322 is associated 
with each set of computation elements 308. Here again, 
the accumulator 320 associated with each set accumulates 
a total of the values of A(k,i,j) over successive cycles 
or successive values of k. Because these totals for 
different sets are accumulated with different values of 
i, they are referred to by the notation SAi. Likewise 
the accumulator 322 associated with each set accumulates 
a total of the values of 2<k,i,j) over successive cycles, 
referred to by the notation SSi. 

[0050] Log circuits 324 and 326, and difference and 
multiplication circuits 328 are also provided. These are 
identical to the corresponding elements 124, 126 and 128 
discussed above with reference to Fig. 3a. Circuits 328 
yield an N 2 element vector 43 of the form 2 0 log ( 2 Ai /E S i ) ... 
201og (SAi/SSi) ... 2 0 1 og ( S A N2 / 2 S N2 ) . A comparator 332 

receives this vector and selects the element N{i) having 
the lowest value. This implicitly selects the value of i 
which yields such minimum value. The comparator outputs 
the selected value I. 
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[0051] Determination of I and J in this manner provides 
information which completely specifies the polarization 
of the incoming signal. The values of I and J can be 
converted directly to (j) and y values corresponding to 
those discussed above with reference to Fig. 1. Thus, 
<j>= J (Ni/3 60 ) and y=I(N 2 /180). These values in turn can be 
converted to the phase angle a, amplitude ratio b/a and 
tilt angle t as discussed above with reference to Fig. 1. 
[0052] The parallel processing system can compute the 
polarization in minimal time. In the particular 

embodiment discussed above, the first group of 
computation elements and associated accumulators provide 
the totals SAj and SSj in a number of clock cycles equal 
to the number of sample values in each series, after 
whatever number of clock cycles are required to initially 
fill the shift registers used in this stage. The same 
number of clock cycles are required for the second stage. 
Added to this are the relatively few clock cycles 
required for operation of the elements used to transform 
the accumulator totals into the vectors 39 and 43, and 
for operation of the comparator. 

[0053] Moreover, operation in the digital domain allows 
for parallel operation with inexpensive components in a 
compact arrangement. The various components of the 

polarimeter 17 can be embodied as one or a plurality of fixed 
gate arrays (FGAs) , field programmable gate arrays (FPGAs), 
application specific integrated circuits (ASICs), digital signal 
processors (DSPs) and the like. 

[0054] Although the various components of the 
polarimeter have been shown and described above as 
separate hardware elements for ease of understanding, 
this is not essential. For example, the same physical 
structures can be used as calculation elements and 
accumulators in both stages of operation. For example, 
the calculation elements can be "soft-connected" or 
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connected through controllable switching devices to the 
buffers, so that the inputs of the calculation elements 
can be reconnected as desired during different stages of 
operation. Also, the calculation elements such as 

adders, as well as the accumulators, can be elements of a 
programmable general-purpose device, and the connections 
required to move the sample values can be made using 
appropriate software instructions for routing the data 
within such a device. 

[0055] The embodiment illustrated above conducts all of 
the various sets of calculations and all of the 
accumulations required in each stage in parallel with one 
another. This can be comprised as, for example, by 
splitting the first stage into two or sub-stages 
performed seriatim. Within each sub-stage, some of the 
calculations and accumulations required in the first 
stage are performed in parallel. The number of clock 
cycles required for the first stage is multiplied by the 
number of sub-stages, but hardware requirements are 
reduced. The same approach can be applied to the second 
stage . 

[0056] Conversely, the stages can be combined with one 
another. A two-dimensional array including NixN 2 sets of 
calculation elements and accumulators can be used in a 
single stage to provide a two-dimensional vector with 
elements corresponding to all possible values of I and J, 
and the element in the vector representing the best null 
(lowest value of 2 Olog ( SA/SS ) ±j ) can be selected. 
[0057] The calculation elements can use transfer 
functions other than those discussed above. Numerous 
analog polarimeters using different transfer functions 
are known in the art. An array of such polarimeters can 
be simulated using digital computations in parallel 
channels, with different parameters in each channel, in 
the same manner as described above. 
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[0058] Although the invention herein has been described with 
reference to particular embodiments, it is to be understood that 
these embodiments are merely illustrative of the principles and 
applications of the present invention. It is therefore to be 
understood that numerous modifications may be made to the 
illustrative embodiments and that other arrangements may be 
devised without departing from the spirit and scope of the 
present invention as defined by the appended claims. 
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